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基于 样 方 数据 的 云南 松林 分 生长 模型 研究 
一 一 以 云南 省 保山 市 杨柳 白族 琶 族 乡 为 例 
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摘要 : 以 云南 省 保山 市 杨柳 白族 莱 族 乡 实地 调查 的 86 块 云南 松 样 地 数据 ， 使 用 非 线性 拟 合 的 方法 拟 合 优 
选 常用 树木 生长 方程 ， 建 立 了 包括 地 位 指数 、 密 度 指 数 、 平 均 直 径 和 蓄积 量 的 生长 模型 。 经 验 生 长 方程 
Schumacher 的 拟 合 度 与 其 它 方程 相近 ， 但 模型 中 参数 变异 系数 均 比 其 余 方程 低 ， 为 地 位 指数 、 平 均 直 径 和 
蓄积 量 的 最 优 方程 。 
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A Study of Stand Growth Model for Pinus yunnanensis 
( Pinaceae) Based on Plots Data 


一 一 人 Case Study in Yangliu Township, Baoshan, Yunnan Province 
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(1 Centre for Mountain Ecosystem Studies, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650204, China; 
2 Graduate University of Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: Pinus yunnanensis is one of most important timber species in Yunnan Province, and widely distributed in 
southwest China. Studies on growth model have been reported, however, most of which focused on a specific part of 
growth model. To build detailed, easily used and accurate empirical stand growth model from the same dataset, a 
case study was conducted in Yangliu Township, Baoshan, Yunnan Province. A total of 86 sample plots data were 
collected using two stages sampling design. Several popular non-linear growth functions were fitted and compared, 
including Chapman-Richards, Mitscherlich, Schumacher, Gompertz, Logistic, Korf and Allometric function. Mod- 
els of site index, density index, average diameter at breast height (DBH) and stock volume growth model were fitted 
respectively. The different models performed more or less similarly in terms of coefficients of determination and root 
mean square error (RMSE). However, empirical growth function “Schumacher” had lower coefficient of variation 
for all parameters than other models. Schumacher function was the most suitable one for site index, average DBH 
and stock volume growth model in all alternative functions. 


Key words: Pinus yunnanensis; Stand growth model; Nonlinear curve fitting 


云南 松 ( Pinus yunnanensis) 广泛 分 布 于 我 南 省 分 布 面积 约 5x10" hm ， 占 云南 省 森林 面积 





国 西南 地 区 ， 是 云南 省 主要 用 材 树种 之 一 。 在 云 





米 


米 米 























基金 项 目 ， 农村 生物 能 源 生产 项 目 (0704711111) 和 “十 一 五 ”国家 科技 支撑 计划 亚热带 森林 区 营 林 固 碳 技术 研究 示范 项 目 





(2008BAD95 B09) 
通讯 作者 : Author for correspondence; E-mail: jxu@ mail. kib. ac. cn 
收 稿 日 期 2010-11-05 ，2011-01-20 接受 发 表 














作者 简介 : 郎 荣 (1984-) 女 ， 硕 士 研究 生 ， 主 要 从 事 森 林 生 长 模型 和 植被 遥感 方面 的 工作 。E-mail: langrong@ mail. kib. ac. cn 


DOI: 10.3724/SP.J.1143.2011.10192 


的 52% ， 鞭 积 量 约 3x10 m?， 为 云南 省 有 林地 
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蓄积 总 量 的 32% ( 雇 声 巾 等 , 2009 ) 。 建 立 林 分 
生长 模型 模拟 和 预测 林 分 的 生长 动态 对 于 森林 的 
经 营 和 管理 具有 重要 意义 。 林 分 生长 模型 是 描述 
林 分 生长 与 林 分 状态 和 立地 条 件 关系 的 一 个 或 一 
组 数学 函数 模型 。 已 有 一 些 研 究 从 不 同方 面 对 云 
南 松 的 生长 动态 进行 了 研究 ， 包 括 林 分 高 度 结构 
( 李 本 德 等 , 1984) 、 林 分 树 高 和 胸径 的 生长 关系 
( 张 学 权 和 喻 朝阳 , 2000) 、 林 分 系统 的 自然 演变 
( 杨 永 祥 等 , 1991)、 自 然 稀 蚊 模型 ( 江 希 铀 等 ， 
2001; 吴 承 祯 等 , 2005 ) 、 径 阶 分 布 模型 〈 李 一 
清 , 1992; 姜 舌 等 , 2008) 、 林 分 经 营 密度 ( 杨 
再 强 和 谢 以 萍 ，1998 ; 卢 学 苏 , 2000) 、 林 分 因子 
的 生长 量 等 ( 雇 声 昭 等 , 2009 ) 。 在 全 林 分 生长 
模型 方面 ， 付 小 勇 (2006) 以 一 平 浪 林场 的 小 班 
数据 系统 地 研究 了 云南 松 的 全 林 分 生长 模型 ， 陈 
绍 安 (2007) 使 用 相同 数据 并 基于 Richard 方程 
研究 了 云南 松 多 形 树 高 生长 模型 。 而 陶 方 红 
(2007) 对 Mitscherlich 生长 模型 进行 了 讨论 。 

全 林 分 生长 模型 使 用 主要 的 林 分 特征 因子 模 
拟 林 分 总 体 或 平均 单 株 木 的 生长 动态 ， 模 型 的 建 
立 和 使 用 较为 简易 ， 特 别 是 对 人 工 林 模 拟 的 通用 
性 和 准确 性 均 较 好 ( Vanclay, 1994)。 系 统 的 云 
南 松 全 林 分 生长 模型 研究 还 很 少 ， 付 小 勇 
(2006) 建立 的 基于 小 班 数据 的 单 形 生长 模型 拟 
合 度 低 ， 尽 管 多 形 生长 模型 拟 合 度 高 ， 但 多 形 地 
位 指数 的 计算 较为 困难 ， 在 实际 应 用 时 存在 诸多 
不 便 。 此 外 林 分 生长 模型 的 研究 以 样 地 数据 特别 
是 固定 样 地 数据 为 主 ， 由 于 实测 样 地 数据 难以 获 
得 ， 建 立 简单 易 用 和 可 靠 的 全 林 分 生长 模型 仍 有 
许多 工作 值得 去 做 。 本 文 以 实地 调查 获得 的 第 一 
手 云南 松 样 方 数据 为 基础 研究 了 云南 松 的 经 验 生 
长 模型 。 立 地 质量 和 林 分 密度 是 影响 林 分 生长 的 
重要 因子 ， 引 入 立地 质量 和 密度 的 指标 可 以 更 好 
的 反映 不 同 立地 条 件 和 林 分 密度 下 的 林 分 生长 ， 
使 模型 具有 更 好 的 通用 性 。 


























1 研究 资料 与 方法 
1.1 研究 区 域 概述 

杨柳 白族 苏 族 乡 位 于 云南 省 保山 市 隆 阳 区 ， 东 经 
98°50'12.6”-99°5'30.6”， 北 纬 25°4'4. 1”-25°20'14. 4”， 
面积 474.6 km 。 地 处 低 续 高原， 海拔 差异 较 大 ， 地 形 地 
貌 多 样 。 隆 阳 区 属于 山地 中 亚热带 季风 气候 ， 年 平均 气 


















































温 16.3%C， 年 降雨 量 1037 mm ( 李 赤 , 2009) 。 依 据 2004 
年 保山 市 杨柳 白族 苏 族 乡 二 类 林业 调查 小 班 数据 ， 农 地 
为 主要 土地 利用 类 型 ， 森 林 和 覆盖 率 为 33% 。 其 中 云南 松 
纯 林 面积 7969 hm ， 占 全 部 森林 面积 的 50.48% ; 以 云南 松 
为 优势 种 的 混交 林 面 积 873 hm*， 为 森林 面积 的 5.30% ， 
现 有 云南 松林 多 为 近 熟 林 ， 无 过 熟 林 。 成 熟 林 在 被 划分 
为 水 土 保持 、 水 源 涵养 或 环境 保护 的 国家 、 地 方 生态 公 
益 林 里 保存 较 好 。 一 般 用 材 林 中 未 被 砍伐 利用 的 成 熟 林 
极 少 ， 中 龄 林 则 多 为 陆续 砍伐 后 自然 更 新 形成 。 林 分 以 
人 工 起 源 为 主 ， 部 分 林 分 为 上 世纪 80 年 代 飞 播 形成 。 
1.2 样 方 调查 

样 方 调查 采用 两 阶段 抽样 的 方法 。 第 一 阶段 抽样 以 
林业 二 类 调查 的 小 班 数据 为 参考 ， 选 取 所 有 云南 松 纯 林 
小 班 以 龄 组 分 层 。 龄 组 依据 2004 年 云南 省 森林 资源 规划 
设计 调查 操作 细则 中 的 标准 划分 ，10 年 及 以 下 林 龄 的 林 
分 为 幼 龄 林 ，11 ~ 20 年 林 分 为 中 龄 林 ，21 ~ 30 年 林 分 为 
近 熟 林 ，31 ~ 50 年 林 分 为 成 熟 林 ，51 年 及 以 上 林 龄 的 为 
过 熟 林 。 每 层 随机 抽取 6 个 小 班 作为 抽样 底 图 。 每 个 小 
班 内 设置 2 个 20mx20m 的 样 方 ， 样 方 间 至 少 间隔 50 m。 
第 一 阶段 抽样 于 2009 年 5-6 月 在 保山 市 林业 局 林业 技 
术 推 广 总 站 协助 下 完成 。 为 了 选取 预定 数目 的 成 熟 林 样 
地 ， 也 在 杨柳 乡 附近 的 瓦 房 乡 调查 了 9 个 样 地 。 根 据 第 
一 阶段 样 方 数据 的 平均 值 和 标准 差 计算 容忍 误差 为 20% 
时 各 龄 组 所 需 样 本 量 ( Krebs, 1999)， 据 此 实行 第 二 阶 
段 抽样 ， 补 充 各 龄 组 所 需 样 方 ， 第 二 阶段 样 方 调查 于 
2010 年 7 月 在 杨柳 白族 苏 族 乡 林业 站 的 协助 下 完成 。 两 
阶段 共 实 地 调查 了 91 个 样 方 ， 包 括 幼 龄 林 5 个 样 方 ， 中 
龄 林 17 个 样 方 ， 近 熟 林 43 个 样 方 ， 成熟 林 26 个 样 方 ， 
无 过 熟 林 样 方 。 其 中 幼 龄 林 和 中 龄 林 样 方 数量 未 满足 两 
阶段 抽样 的 数目 要 求 ， 主 要 由 于 近年 来 造林 较 少 选用 云 
南 松 ， 自 然 更 新 的 林 分 常 林 相 不 齐 ， 林 木 间 林 龄 差别 较 
大 ， 不 宜 作为 同龄 林 生 长 模型 的 样 地 。 总 体 来 说 ， 所 调 
查 样 地 林 分 群落 结构 相对 简单 ， 偶 有 样 地 乔木 层 出 现 华 
山 松 ， 仅 有 少数 生态 公益 林 样 地 林 下 灌木 层 发 达 ， 物 种 
组 成 较为 丰富 ， 其 余 样 地 由 于 放牧 干扰 等 原因 林 下 灌木 
和 草本 稀 琉 。 研 究 区 域 的 云南 松林 及 样 地 分 布 如 图 1 所 
示 ， 其 中 12 个 样 地 由 于 天 气 等 原因 无 GPS 记录 ， 未 能 
够 在 图 上 显示 。 

调查 时 记录 样 地 基本 信息 和 各 林 分 因子 测量 值 。 样 
地 基本 信息 包括 GPS 坐标 、 海 拔 、 坡 度 、 坡 向 和 林 龄 等 
信息 。 云 南 松 每 年 在 树 的 顶端 新 生 轮 生 枝 一 轮 ， 通 过 查 
轮 生 枝 可 获得 林 龄 ， 坡 度 使 用 六 五 式 指南 针 测量 。 测 量 
的 林 分 因子 包括 胸径 和 树 高 ， 胸 径 为 每 木 检 尺 ， 使 用 胸 
径 卷 尺 测量 。 我 国 国家 森林 资源 连续 清查 技术 规定 和 云 
南 省 森林 资源 规划 设计 调查 操作 细则 规定 乔木 检 尺 起 测 
胸径 为 5 cm。 本 研究 里 幼 龄 林 中 大 量 树木 胸径 小 于 起 测 
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胸径 ,为 了 更 准确 的 反映 幼林 的 生长 状况 ， 对 所 有 幼 龄 
林 的 云南 松 胸径 进行 了 测量 。 其 余 龄 组 中 胸径 小 于 5 cm 
的 树木 常 为 被 压 木 或 自然 更 新 的 幼 树 ， 其 蕃 积 量 所 占 比 
例 非 常 小 ， 为 了 与 林业 调查 的 起 测 胸 径 一 致 ， 对 这 些 龄 
组 林 分 的 胸径 测量 仅 针对 5 cm 及 以 上 的 树木 。 树 高 则 根 
据 径 级 分 布 情况 选取 10 ~ 1 株 测 量 ， 包 括 3 株 最 高 的 林 
木 树 高 作为 优势 高 。2009 年 调查 时 使 用 伸缩 式 测 高 杆 测 
量 ，2010 年 调查 时 使 用 超声 波 测 高 仪 “ Forestor Vertex” 


测量 。 
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图 1 杨柳 乡 云南 松林 及 样 地 分 布 图 


Fig.1 Distribution of Pinus yunnanensis stands and location 


of sample plots in Yangliu Township 


1.3 数据 预 处 理 

计算 和 整理 各 样 地 的 平均 直径 、 优 势 高 、 平 均 高 、 
林 分 密度 和 蓄积 量 。 平 均 直 径 为 所 有 检 尺 云南 松 胸径 的 
算术 平均 值 ， 优 势 高 为 样 地 内 3 株 最 高 优势 木 树 高 的 平 
均值 ,平均 高 依据 指数 式 方程 的 数 式 法 代入 平均 直径 计 
算 ， 计算 密度 时 对 于 坡度 大 于 5 度 的 样 地 进行 坡度 校 
正 。 单 株 蓄积 量 依 据 中 华人 民 共 和 国 农林 部 部 颁 标准 
LY208-77 (1978) 中 云南 松 的 二 元 材积 式 计 算 ( 式 1)， 
林 分 蕾 积 量 为 单 株 蓄积 量 和 密度 的 乘积 。 























V=0. 0000582901 17 时 DD! 9796344 。 H" 90715154 ( 1 ) 


检查 所 有 样 地 数据 后 剔除 了 5 个 样 地 数据 ， 其 余 86 
个 样 地 数据 参与 模型 的 拟 合 。 押 除 的 样 地 中 1 个 样 地 经 
确认 为 思茅 松 ，2 个 为 林 相 不 齐 无 法 确定 林 分 的 林 龄 ， 
其 余 2 个 可 能 为 测量 或 记录 错误 ， 其 计算 的 林 分 因子 值 
与 其 余 样 地 数据 偏离 较 远 。 主 要 林 分 因子 的 平均 值 、 最 
大 值 、 最 小 值 和 标准 差 列 于 表 1 。 






































表 1 主要 林 分 因子 的 统计 描述 


Table 1 Descriptive statistics of main stand factors 




















去 首 关 
林 分 因子 平均 值 ”最 大 值 最 小 值 Ri 和 
Stand factors Mean Maximum Minimum ,)... 
division 
林 龄 23.9 41.0 10.0 7.8 
Stand age (year) 
分 辣 诺 
林 分 密度 . 5 1607.8 5585 238 1093.7 
Stand density ( N/hm’ ) 
平均 胸径 
3 
Mean DBH (em) 13.8 30.8 4.3 5.0 
1 18.1 4.0 3.2 
Mean tree height (my) 了 8 
人 13.3 20.5 5.3 3.5 


Dominant tree height (my) 


1.4 备 选 生长 方程 

本 文选 用 经 验 生长 方程 Schumacher， 理 论 生长 方程 
Chapman-Richards 、 Mitscherlich、 Gompertz、 Logistic 、 Korf 
和 Allometric 作为 生长 模型 的 备 选 方程 ， 方 程 形式 根据 
重 宪 宇 (2006) 和 Origin 内 置 函数 形式 确定 〈 表 2) 。 理 
论 生 长 方程 Chapman-Richards 、Mitscherlich 、Gompertz 和 
Korf 方程 中 的 参数 a 为 树木 生长 的 最 大 值 参数 ，b 为 生 
长 速率 参数 ，Logistic 方程 中 a 为 最 大 值 参数 ，e 为 生长 
速率 参数 。 模 型 的 评价 指标 包括 决定 系数 (R ) 、 均 方 
根 误差 (RMSE) 、 拟 合 参数 的 变异 系数 ( C. V)。 均 方 
根 误差 作为 模型 误差 评价 的 一 个 常用 指标 ， 反 映 观测 值 
和 模型 估 测 值 的 差异 。 各 方程 的 拟 合 使 用 Origin 8 Pro 的 
非 线性 拟 合 工具 完成 。 该 软件 包含 大 量 内 置 非 线性 方程 
可 供 拟 合 ， 这 些 内 置 方程 不 需 给 定 拟 合 参数 的 初始 值 ， 
可 避免 因 人 为 给 定 初始 参数 值 不 适 造成 方程 拟 合 不 收敛 
的 问题 。Origin 8 Pro 的 内 置 非 线性 方程 中 无 Schumacher 
和 Korf 的 方程 ， 因 此 通过 自 建 这 两 个 方程 以 供 拟 合 。 


















































表 2 生长 模型 备 选 方程 列表 
Table 2 Alternative fitting functions for growth model 


备 选 方程 名 称 具体 方程 


Alternative functions 














Name of alternative functions 








(1) Chapman-Richards Y=a (1-e”)° 
(2) Mitscherlich Y=a (1-e™™) 
(3) Schumacher Ye 

(4) Gompertz Wee 
(5) Logistic pe 
(6) Korf Y=ae 

(7) Allometric Y=att 
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2 林 分 生长 模型 
2.1 优势 高 生长 模型 和 地 位 指数 

林 分 树 高 是 同龄 林立 地 质量 指标 中 最 广 为 接 
受 和 使 用 的 指标 (Skovsgaard 和 Vanclay, 2008 ) 。 
地 位 指数 (site index) 指 某 树种 在 基准 林 龄 时 的 
林 分 优势 高 或 平均 高 ， 本 文选 用 优势 高 和 林 龄 计 
算 的 地 位 指数 作为 立地 质量 的 评价 指标 。 以 各 样 
地 的 优势 高 (HH ) 为 因 变量 ， 林 龄 (1) 为 自 变 
量 拟 合 备 选 生长 方程 ， 拟 合 结果 和 评价 指标 列 于 
表 3， 所 有 拟 合 方 程 的 下 检验 均 为 极 显 著 (P< 
0.01)。 

各 备 选 模 型 的 决定 系数 和 均 方 根 误差 非常 接 
近 ， 但 拟 合 参数 的 变异 系数 及 上 检验 显著 性 差异 较 
大 。 各 模型 的 均 方 根 误差 为 2.03 ~2.04， 其 中 以 
Schumacher 方程 的 参数 变异 系数 最 小 ，a 和 b 的 变 
异 系数 分 别 为 5.65% 和 9.36% ， 且 + 检验 均 达到 极 
显著 水 平 ， 综合 考虑 各 评价 指标 后 选取 Schumacher 
为 优势 高 生长 模型 ( 式 2) ， 结 果 如 图 2 所 示 。 
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林 龄 Stand age (year) 
图 2 ”Schumacher 优势 高 模型 


Fig.2 Stand dominant height model based on Schumacher equation 


15. 62437 


H,=25.39585. e (2) 





根据 地 位 指数 (SI) 的 计算 方法 ， 以 20 年 为 基 
准 林 龄 ， 各 样 地 的 地 位 指数 可 由 优势 高 和 林 龄 计 
算得 出 〈 式 3) 。 


15. 62437 15.62437 
# “0 


SI=H,* e' ” (3) 
2.2 密度 指数 

林 分 密度 是 森林 经 营 者 可 调节 和 控制 林 分 生 
长 和 发 育 的 主要 因子 ,将 其 引入 方程 可 预 估 不 同 
密度 林 分 的 生长 状态 ， 使 模型 具有 更 好 的 通用 
性 。 本 文选 用 林 分 密度 指数 ( stand density in- 
dex，SDI) 为 林 分 密度 指标 。Reineke 发 现 完满 
立木 度 、 未 经 间伐 的 同龄 林 中 ， 单 位 面积 株数 
(N) 与 林 分 平均 胸径 〈D) 之 间 呈 寡 函 数 关系 ， 表 
示 为 N=a D-4 (Reineke，1933 ) ， 对 数 转换 后 成 
为 线性 方程 ( 式 4) 。 株 数 密度 和 平均 直径 经 对 
数 转换 后 拟 合 该 式 ， 求 出 方程 中 的 两 个 参数 值 。 


lgN=lga-8lgD (4) 


完满 立木 度 方程 的 决定 系数 RE =0.63，RMSE 
=0.17，lg a=4.64391, B=-1.37327， 两 个 参 
数 的 1 检验 均 达 极 显著 水 平 。 根 据 密度 指数 
(SDI) 的 计算 方法 ,依据 林 分 密度 指数 方程 
( 式 5) 计算 各 样 地 的 密度 指数 ， 本 文 基准 直径 
(D,) 设 定 为 15 cm。 

















D 
SDI=N (一 ) -3 (5) 
D 


表 3 优势 高 模型 拟 合 结果 


Table 3 ” Nonlinear curve fitting results of dominant tree height 








备 选 方程 评价 指标 Evaluation indicators 参数 值 (C.V, %) Parameters 

Alternative functions R2 RMSE a b C 
Chapman-Richards 0.68 2.03 22.36816™ (20.79) 0.04301 (58.17) 1.22543 ™ (36.52) 
Mitscherlich 0.68 2.03 25.6737™ (13.52) 0.02931 (21.02) 

Schumacher 0.67 2.03 25.39585™ (5.65) 15.62437™ (9.36) 

Gompertz 0.68 2.03 20.10305 (10.11) 0.07049 (24.67) 12.11452™ (11.52) 
Logistic 0.68 2.03 18.83803 (7.19) 5.42395 ™ (27.08) 0.10416™ (19.24) 
Korf 0.68 2.03 44.60177 (80.10) 7.08325 * (48. 23) 0.55093 (64.39) 
Allometric 0.67 2.04 1.47492 (19.49) 0.68™ (8.57) 


“为 P<0.05, ”为 P<0.01 
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2.3” 林 分 平均 直径 生长 模型 

本 研究 中 林 分 平均 直径 随 林 龄 增长 呈 近 似 线 
性 的 增长 。 以 平均 直径 为 因 变 量 ， 林 龄 为 自 变量 
拟 合 备 选 生长 方程 ， 拟 合 结果 如 表 4 所 示 ， 各 备 
选 方程 下 检验 均 达 极 显著 性 水 平 (P<0.01) 。 

各 模型 的 决定 系数 和 均 方 根 误差 非常 接近 ， 
决定 系数 在 0.57 至 0.61 之 间 ，RMSE 在 3.20 ~ 
3.34 之 间 ， 模 型 估 测 值 和 观测 值 差 异 较 小 。Korf 
方程 的 决定 系数 最 高 ，Schumacher 方程 的 决定 系 
数 最 低 。 比 较 各 模型 的 参数 变异 系数 和 显著 性 水 
平 ， 则 以 Schumacher 方程 表现 最 好 。 理 论 牛 长 
方程 Chapma-Richards 、Mitscherlich 、Gompertz 和 
Korf 方程 的 参数 a 为 林木 生长 的 最 大 值 参 数 ， 但 
拟 合 结果 中 的 参数 值 已 不 再 符合 实际 的 平均 直径 
值 ， 使 得 拟 合 的 理论 生长 方程 失去 了 其 生物 学 意 
义 ; 同时 这 几 个 方程 的 参数 变异 系数 也 很 大 。 
Schumacher 方程 为 平均 直径 生长 模型 的 最 优 方程 
( 式 6) ， 拟 合 曲线 如 图 3 所 示 。 
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图 3 Schumacher 平均 直径 生长 模型 


Fig.3 Average DBH growth model based on Schumacher equation 























22. 31129 





D=33.92571. e- (6) 


林 龄 单 因 变量 的 平均 直径 生长 模型 的 拟 合 度 
较 低 ， 引 入 地 位 指数 和 密度 指数 再 次 参数 化 生长 
模型 。 最 大 值 参 数 a 表达 为 地 位 指数 的 函数 ， 生 
长 速率 参数 表达 为 密度 指数 的 函数 ， 均 采用 窜 函 
数 的 形式 。 以 各 样 地 的 平均 直径 为 因 变 量 ， 林 
龄 、 地 位 指数 和 密度 指数 为 自 变量 拟 合 再 次 参数 
化 的 Schumacher 方程 ( 式 7)。 


a,* (SDI/1000)” 


y=a1* (SI)”".e-— 





(7) 


非 线 性 拟 合 的 结果 为 R* =0.69，RMSE =2. 84， 
al =5.75854, a,=23.09328, b,=0.72941, b,= 
0.21198， 拟 合 参 数 t 检验 中 a 达 显著 性 水 平 ， 
其 余 均 达 极 显著 水 平 。 林 龄 单 因 变量 的 平均 直径 
生长 模型 引入 地 位 指数 和 密度 指数 后 ， 模 型 的 拟 
合 度 略 有 提高 。 模 型 的 决定 系数 由 0. 57 提高 至 
0. 69 ， 均 方 根 误差 由 3.34 降 至 2. 84， 新 变量 的 
引入 提高 了 模型 的 解释 度 但 提高 幅度 不 大 。 包 括 
林 龄 、 地 位 指数 和 密度 指数 的 平均 直径 生长 模型 
如 式 8 所 示 。 
































”23.09328 ， (SDH1A1000) 2 


D=5.75854 . (SD)°™™.e 1 





(8) 
2.4 林 分 蕃 积 量 生 长 模型 
以 林 分 蓄积 量 为 因 变 量 ， 林 龄 为 自 变 量 ， 拟 
合 备 选 生长 方程 的 结果 列 于 表 5， 所 有 方程 的 下 
检验 均 为 极 显著 (P<0. 01)。 














表 4 ” 林 龄 单 因 变量 的 林 分 平均 直径 生长 模型 拟 合 结果 


Table 4 Fitting results of stand average DBH growth model with age variable 








备 选 方 程 评价 指标 Evaluation indicators 参数 值 (C.V，% ) Parameters 

Alternative functions R2 RMSE b 2 
Chapman-Richards 0.60 3.22 45454. 15001 (502395) 0.000007 (531765) 0.94554 (53.35) 
Mitscherlich 0.60 3.20 —125617. 54884 (163101) -0.000004 (163127) 

Schumacher 0.57 3.34 33.92571 ™ (9.77) 22.31129™ (11.87) 

Gompertz 0.60 3.21 144.48037 (299. 80) 81.4220 (189.74) 0.01557 (131.34) 
Logistic 0.60 3.22 70.90022 (184.14) 13. 8$122 (171.03) 0.04516 (50.93) 
Korf 0.61 3.20 1.57904 (244.19) -0.51212 (278.42) -0.44083 (114.77) 
Allometric 0.60 3.20 0.63603 (31.46) 0.94546 (9. 84) 


“为 P<0.05, ”为 P<0.01 
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表 5 鞭 积 量 与 林 龄 单 因 变量 生长 模型 拟 合 结果 

Table 5 Fitting results of the stand volume growth models with age variable 
备 选 方程 评价 指标 Evaluation indicators 参数 值 (C.V, % ) Parameters 
Alternative functions R? RMSE a b C 
Chapman-Richards 0.37 64.05 41194. 28514 (18340.41) 0.00050 (14224.09) 1.3398 (102.66) 
Mitscherlich 0.37 63.62 —158.13209 (78.14) -0.02179 (56.40) 
Schumacher 0.36 64.11 487.32016 “(22.76) 34.85857 ™ (18.86) 
Gompertz 0.37 64. 18 901.62656 (266. 81) 0.02554 (143.46) 53.48375 (171.77) 
Logistic 0.37 64.29 852.61786 (354.34) 27.9161 (306.38) 0.05756 (75.02) 
Korf 0.37 64. 04 0. 00363 (7669.97 ) -6.88359 ( -1061.55) -0.12718 (725. 42) 
Allometric 0.37 64.04 1.59752 (74.57) 1.33071 ™ (16.36) 


“为 P<0.05, ”为 P<0.01 











林 龄 单 因 变量 的 蓄积 量 生 长 模型 决定 系数 
低 ， 均 方 根 误差 值 高 ， 各 模型 的 估 测 值 与 观测 值 
间 存 在 差异 大 ， 林 龄 单 因 变量 对 蓄积 量 生 长 的 解 
释 程 度 有 限 。Schumacher 模型 的 决定 系数 略 低 ， 
但 拟 合 参 数 变异 系数 最 低 ， 参 数 a 和 b 的 变异 系 
数 分 别 为 22.76% 和 18. 86% ， 检验 均 达 极 显著 
水 平 ， 因 此 选择 Schumacher 方程 为 蔓 积 量 生 长 
的 最 优 模型 ( 式 8) 。 从 各 备 选 方程 的 拟 合 参数 
来 看 ，Chapman-Richards 、Mitscherlich 、Gomper- 
tz、Logistic 和 Korf 模型 的 所 有 拟 合 参 数 的 1 检验 
均 为 不 显著 ，Allometric 方程 的 一 个 拟 合 参数 也 
未 达 显 著 性 水 平 。 理 论 生 长 方程 如 Chapman- 
Richards 、Mitscherlich 的 最 大 值 参 数值 已 失去 应 
有 的 生物 学 意义 。 
























































34. 85857 


SV=487.3202.e6 





(9) 


同 理 引 入 地 位 指数 和 密度 指数 再 次 参数 化 
Schumacher 模型 ( 式 7)， 非 线性 拟 合 后 的 蓄积 量 
生长 模型 R =0. 89，RMSE = 26. 39，ai =41. 76061 ， 
au =45.44669，b| =1.08456，b, = -0.81342， 所 
有 拟 合 参数 的 1 检验 均 达 显著 性 水 平 。 林 龄 单 因 
变量 的 蓄积 量 生长 模型 引入 地 位 指数 和 密度 指 
数 后 ,， 模 型 的 拟 合 度 得 到 明显 提高 ,决定 系 数 
由 0.36 提高 至 0.89， 均 方 根 误差 由 64. 11 降 至 
26.39， 引 入 的 新 变量 对 模型 解释 度 的 提高 贡献 
大 。 包 含 林 龄 、 地 位 指数 和 密度 指数 的 Schuma- 
cher 林 分 蓄积 量 模型 如 式 10 所 示 。 


_45.44669 + ( SD//1000) -0.81342 
t 


























3 结论 与 讨论 

林 分 优势 高 、 平 均 直 径 和 蓄积 量 的 备 选 生长 
模型 拟 合 时 ， 各 方程 的 决定 系数 和 均 方 根 误差 比 
较 接近 ,但 拟 合 参数 的 变异 系数 及 1 检验 结果 差 
别 较 大 ， 均 以 Schumacher 的 变异 系数 最 小 ， 建 
立 的 模型 稳定 性 好 。 本 研究 中 Schumacher 方程 
的 通用 性 较 其 余生 长 方程 更 好 ， 这 与 付 小 勇 
(2006) 研究 中 得 出 的 结论 一 致 。 理 论 生长 方程 
虽 能 够 对 方程 中 的 参数 给 出 合理 的 生物 学 解释 ， 
但 在 本 研究 中 该 类 方程 拟 合 效果 不 佳 ， 直 接 使 用 
某 一 生长 方程 建立 的 生长 模型 可 能 并 不 是 最 合适 
的 模型 。 引 入 地 位 指数 和 密度 指数 再 次 参数 化 模 
型 能 够 大 幅 提 高 蓄积 量 生 长 模型 的 拟 合 度 ， 但 对 
平均 直径 生长 模型 的 提高 幅度 较 小 ， 主 要 原因 为 
地 位 指数 和 密度 指数 与 蓄积 量 、 平 均 直 径 的 相关 
性 程度 不 同 。 

各 模型 的 拟 合 优选 过 程 中 通过 均 方 根 误差 指 
标 对 模型 的 估计 误差 进行 了 评价 ， 林 分 优势 高 的 
估计 值 与 实测 值 差异 为 2.03， 包 含 地 位 指数 和 
密度 指数 的 平均 直径 和 蓄积 量 生长 模型 的 差异 分 
别 为 2. 84 和 26.39。 蓄 积 量 生长 模型 的 均 方 根 
误差 略 高 ， 其 生长 可 能 受到 除 林 龄 、 立 地 和 密度 
外 其 它 因素 的 影响 。 本 研究 使 用 同 组 数据 对 云南 
松林 分 的 地 位 指数 、 密 度 指数 、 平 均 直 径 和 蓄积 
量 模型 进行 了 研究 ， 建 立 的 模型 拟 合 度 较 基于 小 
班 数据 的 单 形 生长 模型 ( 付 小 勇 , 2006) 更 好 ; 
与 云南 松 的 多 形 生 长 模型 ( 陈绍 安 , 2007) 相 
比 ， 本 研究 中 模型 的 拟 合 度 略 低 ， 但 模型 形式 简 
单 ， 便 于 实际 应 用 。 






































SV=41. 76061 + (SI) "6 .e 
(10) 


本 研究 对 模型 的 估计 误差 进行 了 评价 ， 其 准 
确 性 可 满足 研究 区 域内 云南 松林 分 生长 模拟 的 需 
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要 ， 适 用 于 模拟 云南 松 同龄 纯 林 林 分 的 生长 。 由 
于 大 量 野外 样 方 调查 受 时 间 和 经 费 限 制 ， 未 使 用 
独立 的 数据 进行 验证 ， 模 型 在 其 它 区 域 的 适用 性 
尚 有 待 进一步 验证 。 研 究 区 域内 无 云南 松 过 熟 
林 ， 因 此 生长 模型 的 拟 合 数据 中 缺乏 此 类 样 地 数 
据 ， 故 对 50 年 或 更 长 时 间 的 林 分 预测 可 能 与 实 
际 情况 存在 差别 。 雇 声 周 等 (2009) 对 云南 松 
生长 量 的 研究 得 出 25 ~ 50 年 期 间 直 径 生 长 旺盛 ， 
处 于 材积 速生 阶段 ，50 年 以 后 达到 速生 后 期 ， 
各 项 生长 量 开始 出 现 明显 下 降 ， 因 此 补充 50 年 
左右 和 更 大 林 龄 的 样 地 数据 对 于 确定 林 分 生长 量 
开始 下 降 的 拐点 非常 重要 。 
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